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基于 小 波 阔 值 和 全 变 分 模型 的 图 像 去 噪 - 
ik oh, ARRI 


(西安 邮电 大 学 自动 化 学 院 , 8710121) 


的 去 噪 方法 。 利 用 小 波 变换 的 时 频 域 特性 将 含 唆 图 像 分 解 得 到 各 维度 小 波 系数 ,对 低频 小 波 系数 利用 全 变 分 模型 去 噪 ， 
对 于 高 频 系数 根据 不 同 分 解 尺度 选择 不 同 的 最 佳 阅 值 去 噪 ， 克 服 了 统一 赋值 的 不 足 ， 增 强 了 碍 法 的 自 适 应 性 。 理 论 分 
析 和 仿真 实验 结果 表明 ， 所 提 方 法 兼顾 了 小 波 变 换 和 全 变 分 模型 的 去 噪 优点 ， 在 有 效 去 除 噪声 的 同时 更 完整 地 保留 了 
图 像 的 边缘 和 细节 信息 ， 有 较 高 的 结构 相似 度 和 烽 值 信 骂 比 。 
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Method for image denoising based on wavelet transform and total variational model 


Zhang Hong, Zhou Xiaolif 


(School of Automation, Xi’ an University of Posts & Telecommunications, Xi'an 710121, China) 


Abstract: As the drawbacks of traditional wavelet threshold algorithm for image noise reduction that it could not eliminate the 
noise thoroughly and blur the image edge, this paper proposed an image denoising method based on wavelet transform and total 
variation model. Firstly, the wavelet coefficients of each dimension separated from the noisy image through the wavelet 
transform, then use the total variation model to filter the low frequency coefficients, and for the high frequency coefficients, 
different thresholds in different decomposition scales for eliminating the noise overcame the lack of the universal threshold and 
enhanced the adaptability of the algorithm. Finally, the wavelet inversion obtained the denoised image. Theoretical analysis 
and Experimental results show that this algorithm with the advantages of both wavelet transform and total variation model, it 
can remove noise effectively, keep image edges and details very well, the structural similarity and peak signal-to-noise ratio 
have also improved. 
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差 较 大 。 为 解决 该 问题 , HE T Ae EMTS, 
但 在 保留 图 像 细节 时 仍 有 一 定 缺 陷 , 如 基于 邻 域 小 波 系数 的 图 
图 像 在 获取 、 传 输 、 转 换 等 过 程 中 不 可 避免 会 受到 噪声 的 像 去 噪 方法 中， 虽然 研究 了 领域 小 波 系 数 对 闵 值 选取 的 影响 ， 

干扰 ， 对 图 像 的 识别 、 分 析 和 理解 带 来 不 利 影响 册 。 因 此 需要 改进 了 统一 阔 值 ， 但 没有 考虑 子 带 对 阔 值 的 影响 。 利 用 拓 模 滤 
对 图 像 进行 去 噪 处 理 ， 使 后 续 的 分 析 和 信息 提取 更 准确 。 波 器 对 小 波 分 解 系数 进行 自 适 应 滤波 的 方法 09, 但 在 优化 过 程 
前 , 图 像 去 噪 算法 主要 有 :基于 滤波 的 去 噪 算法 、 基 于 多 中 计算 复杂 度 较 高 。 根 据 加 入 阔 值 调节 因子 ， 利 用 采样 点 长 度 
尺度 变换 的 去 噪 算法 和 基于 偏 微分 方程 的 去 噪 算法 证 。 基 于 滤 的 不 同调 节 阔 值 的 方法 00， 由 于 对 阔 值 估计 不 准确 以 及 受 图 像 
波 的 去 噪 算法 包括 中 值 滤波 国 ， 卡 尔 曼 滤波 外 ， 维 纳 滤波 辐 等 ， 大 小 的 影响 ， 去 噪 后 图 像 存 在 边缘 模糊 细节 丢失 的 问题 。 基 于 
这 些 方法 主要 是 通过 滤 除 图 像 的 高 频 部 分 达到 去 噪 的 目的 ， 然 偏 微 分 方程 (partial differential equation,PDE) 的 图 像 去 噪 方法 因 
而 在 滤 掉 噪声 的 同时 会 造成 图 像 的 结构 信息 的 丢失 。 以 小 波 变 其 较 好 的 保 边 性 而 得 到 广泛 应 用 0 区 。Rudin 等 提出 基于 图 像 
换 为 代表 的 多 尺度 变换 具有 和 良好 的 时 频 分 析 特 性 名， 广泛 应 用 的 全 变 分 (total variation,TV)03 模 型 ， 该 方法 是 将 图 像 去 噪 建 模 
于 图 像 去 噪 。 其 中 小 波 阔 值 去 噪 是 小 波 去 噪 的 主要 方法 5 ， 但 成 一 个 能 量 函 数 的 最 小 化 问题 ， 从 全 新 的 角度 解释 图 像 去 噪 的 
该 类 方法 存在 小 波 系数 过 扼杀 现象 ， 导 致 小 波 系数 重 构 时 的 误 过 程 ， 能 够 在 去 噪 的 同时 有 效 保护 图 像 的 边缘 纹理 信息 ， 但 在 
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噪声 较 大 的 情况 下 容易 产 阶梯 效应 05， 从 而 影响 去 噪 效 果 。 

针对 以 上 不 足 ， 利 用 小 波 去 噪 模型 的 时 频 分 析 特 性 和 全 变 
分 模型 保持 图 像 边缘 纹理 等 细节 特征 的 优点 ， 本 文 提出 一 种 自 
适应 小 波 阔 值 和 全 变 分 模型 相 结合 的 图 像 去 噪 方法 。 该 方法 首 
先 对 含 噪 图 像 进行 小 波 变换 ， 将 图 像 分 为 低频 、 高 频 区 域 分 别 
进行 处 理 。 在 图 像 的 低频 区 域 采用 全 变 分 模型 , 对 于 高 频 区 域 ， 
根据 噪声 程度 的 不 同 对 去 噪 力度 进行 自 适应 调整 ， 也 即 采 用 
种 阔 值 随 着 分 解 尺 度 和 子 带 变 化 的 自 适 应 小 波 阔 值 进行 去 噪 ， 
克服 了 统一 阔 值 的 不 足 ， 提 高 了 算法 的 去 噪 性 能 。 此 外 ， 通 过 
对 具有 不 同 噪声 程度 和 不 同 纹理 分 布 的 图 像 进 行 仿真 实验 验证 
本 文 提 出 方法 的 有 效 性 。 


1 ”全 变 分 模型 和 小 波 阅 值 去 噪 


1.1 全 变 分 模型 
全 变 分 模型 首次 将 整体 变 分 引入 到 图 像 去 噪 领域 031。 基 于 
全 变 分 模型 的 图 像 去 噪 是 将 图 像 能 量 用 一 个 函数 表示 ， 在 约束 
条 件 下 ， 通 过 减 小 图 像 系 统 的 能 量 去 除 图 像 噪声 。 因 此 ， 图 像 
的 去 噪 问题 转变 成 了 一 个 能 量 最 小 化 问题 。 
设 为 含 噪 图像 ，Q 表示 图 像 的 定义 域 ， 
foxy) = f(x, y)+nlx, y), HP f(x, y) 为 原始 的 清晰 图 像 ， 
n(x, y) 是 方差 为 o? 的 高 斯 噪声 。 则 f 的 整体 变 分 为 


TVA) = | MW Fad 


图 像 被 噪声 污染 是 一 个 正 问题 ， 根 据 观 测 的 图 像 数 据 、 所 给 的 
限制 条 件 和 所 作 的 假设 来 获得 原 图 像 数 据 。 去 噪 过 程 是 一 个 逆 
问题 ， 可 能 导致 解 的 不 唯一 或 不 稳定 ， 在 数学 上 解决 这 类 问题 
正则 化 方法 对 解 空间 进行 约束 。 因 此 基于 全 变 分 的 图 
果 可 以 归结 为 如 下 最 小 化 问题 : 


l 


(1) 


min E=] Sf) a | AVF Fady © 


项 ， 它 代表 恢复 图 像 和 


其 中 : Sd fady 为 数据 保 


含 噪 图 像 之 间 的 逼近 程度 ， 主 要 起 保留 原 图 像 特性 和 降低 图 像 
le [VIVE P dxdy 为 正则 化 项 


映 出 图 像 中 的 边缘 信息 。 其 导出 的 欧 拉 - 拉 格 朗 日 方程 为 


vi 
| Vf | 


失真 度 的 作 / 它 能 够 很 好 的 反 


Jat fa) =0 (3) 


a 
Ivf | ° 
Vf | 较 大 ， 扩 散 系 数 较 小 ， 因 此 在 边缘 区 扩散 较 弱 ， 从 而 保 
留 了 边缘 。 在 平滑 区 域 ， 扩 散 系 数 较 大 ， 因 此 在 图 像 平 滑 区 域 
的 扩散 能 力 较 强 ， 从 而 去 除 噪声 。 
1.2 i \ RAR 
从 信号 和 噪声 在 小 波 域 中 的 不 同 奇异 性 可 知 09， 有 用 的 信 
号 表现 为 低频 信号 , 高 频 信 号 则 为 噪声 信号 。 经 过 小 波 变 换 后 


从 方程 (3) 可 以 看 出 ， 扩 散 系数 为 在 图 像 边缘 处 ， 


C 
: BEDERA R 


EYONI: R 


fei ARA ERE RRAK, AEAEE 
值 ， 将 绝对 值 小 于 阔 值 的 系数 置 为 零 ， 而 保留 或 收缩 绝对 值 大 


TA 


行 信 
如 图 


波 系 
系数 


时 ， 


法 存 


重 构 


软 阅 什 方法 得 到 的 y, 虽然 连续 怕 


FAW, j 总 是 存在 恒定 的 偏差 


iwi 
2 


$ 


= 


的 


tA 
HA 


区 域 


ANS 


值 的 系数 ， 得 到 估计 小 波 系数 ， 然 后 利用 估计 小 ; 
号 就 是 去 噪 后 的 信号 。 


号 重 构 ， 得 到 的 估计 信 
1 所 示 。 


> 低频 信息 


+ i 


图 1 


小 波 闵 值 去 噪 通 常 又 分 为 硬 闷 值 去 噪 和 软 病 值 去 噪 ， 
硬 病 值 是 将 小 波 系 数 的 绝对 值 不 


数 置 为 0， 
AA Wij MERER 


入 


Wi,j = 


系数 进 
其 去 噪 过 程 


一 > 低 通 滤波 一 去 噪 图 像 


频 信息 


小 波 去 噪 过 程 


W; 如 下 : 


Wj 
0 


>T 


Wij 


<T 


Wij 


其 中 ， 


HBR (ELE AT HE, ANP BELA 
AF BULA DU RBA. KRENT, M 


(4) 


软 阔 值 与 硬 阔 值 不 同 之 处 是 当 小 波 系数 大 于 或 等 于 阔 值 


小 波 系 数值 变 为 该 点 值 与 闵 值 的 差 ， 具体 如 下 : 
A sign(w, Mw; -T) Iw, >T 
Wij = (5) 
0 [w <T 


这 两 种 方法 易于 实现 ， 可 快速 得 到 估计 小 波 系数 ， 但 该 算 


EDE, EEREN 


信号 会 产生 震荡 ， 并 且 会 对 小 波 系数 产生 过 扼杀 现 像 。 


程度 。 


A 


P w ET 处 是 不 连续 的 ， 得 到 的 


Et 好， 但 当 [w |= 时 ， Ww 


了， 影响 了 重 构 信 号 与 真实 信号 的 


自 适应 小 波 阐 值 和 全 变 分 结合 的 去 噪 方法 


针对 统一 阔 值 函数 的 “ 


过 扼杀 ”现象 以 及 造成 图 像 边缘 模 


问题 ， 本 文 提出 一 种 基于 自 适 应 小 波 阔 值 和 全 变 分 模型 相 
波 变换 将 含 噪 图 像 分 为 
。 在 图 像 的 低频 区 域 采 用 全 变 分 模型 去 噪 ， 对 高 频 区 域 采 


的 去 噪 方法 


， 即 通过 小 


2.1 


自 适应 小 波 阐 值 


在 较 小 的 邻 域内 ， 小 波 系数 具有 相关 性 ， 


用 自 适 应 小 波 闵 值 去 噪 ， 最 后 通过 小 波 逆 变换 得 到 去 噪 后 的 医 


氏 频 和 高 频 


因而 在 对 小 波 系 


数 进 


行 闵 值 处 理 时 ， 应 考虑 


像 细 


节 。 因 此 ， 在 闵 值 函数 


其 邻 域 系 数 的 情况 ， 


中 引入 相 邻 系数 07。 


声 信号 的 小 波 系数 集 ， 阔 值 函数 人， 为 


尽 可 能 保留 图 


Hit d, 是 品 
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录用 定稿 
LV 
^ din l-ax— S zaxa 
diry Si (6) 
0 其 他 
1 m=R n=R 3 
其 中 : Se= h Vid), 是 对 各 层 高 频 子 带 的 能 量 估 计 ， 


2 
R m=-R n=-R 


R 是 所 选择 的 邻 域 窗口 大 小 ，m 、nn FIAT A BEE 


的 行 数 和 列 数 ，d,, , 为 为 该 子 带 的 小 波 系 数 , 4 = oy2logN ， 


N 是 信号 的 长 度 ，5 为 噪声 标准 差 。 

小 波 系数 随 分 解 尺 度 的 增加 ， 其 自身 呈现 近似 指数 衰减 ， 
也 即 噪声 在 信号 分 解 的 过 程 中 ， 能 量 逐 渐 减 弱 , 且 分 解 后 的 水 
平和 重 
带 大 ， 说 明 水 平和 垂直 高 频 子 带 中 所 含 噪声 少 。 当 采用 统一 闪 
值 去 噪 时 ， 由 对 角子 带 估计 不 同方 向 的 高 频 子 带 上 的 阐 值 ， 这 
样 会 去 除 掉 垂 直 和 水 平子 带 中 的 部 分 图 像 信息 ， 造 成 图 像 的 模 
糊 。 根 据 信号 和 噪声 在 分 解 中 的 这 一 特性 。 本 文 提 出 通过 变化 
闷 值 函数 式 (6) 中 的 参数 & 来 调整 阔 值 。 根据 小 波 系数 在 图 像 分 


解 中 近似 指数 的 豪 城 特性 ， 取 a= ，; 为 分 解 导数 ， 通 过 调 
整 参数 w ， 可 以 选择 不 同 尺度 上 的 小 波 系 数 对 应 的 阔 值 函数 ， 


这 样 能 够 保证 q ,的 整体 连续 性 ， 从 而 避免 信号 产生 震荡 ， 提 


高 去 噪 效果 。 
2.2， 自 适应 小 波 阔 值 和 全 变 分 相 结合 的 去 噪 方法 

于 小 波 阔 值 函数 去 噪 会 造成 图 像 边 缘 模 糊 ， 而 全 变 分 能 
够 保护 图 像 边缘 , 因此 , 可 以 将 小 波 阔 值 函数 和 全 变 分 相 结合 ， 


7 


达到 去 除 噪声 的 同时 保留 图 像 边缘 的 目的 。 假设 {WV} 是 一 张 


i 


积 空间 ， 构 成 六 (>) 上 的 一 个 多 尺度 分 析 ，8(x) 和 (x) 分 别 


是 相应 的 尺度 函数 和 小 波 函 数 ， 根 据 Mallat 金字 塔 分 解法 ， 图 


像 (x,y) 二 维 离散 小 波 变 换 的 快速 分 解 可 表示 为 


f(xy) =} c, (L(x)+ 


È Lael (LH (3) oe 
k=1,2,3 m,n 
Cmn (LL)=> > g(L WE es, 2m-L 
eZ l,eZ 
a. (LH)=>. > 8 (4 \a(L Verni 2m-l, 
leZbheZ 
A) Ta) P 
leZ heZ 
di (HH)=2 2 a(b h(l) es 2m 
leZ heZ 


其 中 i=1,2,.…,7 为 小 波 分 解 的 尺度 ， 及 (.) , g(.) 分 别 为 高 通 


a a GRAVARE 


VE BE EAR ERR Can (LL) 为 原 图 像 的 近似 表示 ， 


di! 


m,n 


个 方向 的 高 频 小 波 系数 。 


万 的 近似 表示 ， 将 低频 子 带 图 ch, 拓展 至 与 原始 图 像 有 相同 


大 小 ， 这 样 全 变 分 的 保 


Jf chn) dxdy 起 平衡 作 ) 


(LH), dy, (HL), dy, (BEE) 分 别 为 垂直 、 水 平 、 对 角 三 


A IRIT E cv (LL) 是 原 图 像 


a 


I (f -ci n) dxdy 


5 项 就 变 成 


(9) 


其 对 应 的 能 量 泛 函 如 


F: 


A i 2 
Euev= 了 -ceo+| MW Pdxdy (10) 


相应 的 欧 拉 - 拉 格 朗 日 方程 如 下 : 
. Vf 
div( )+A(c,, f) =0 (11) 
| Vf | i 


参数 4 对 扩散 


越 明 显 ， 去 噪 后 的 图 像 就 越 接近 原 图 像 。 本 文 根 据 不 同 图 像 去 


项 [VIVE faddy 和 保 真 项 


1, 4A 值 越 大 时 ， 保 真 项 作用 就 


噪 后 峰值 信 噪 比 的 值 来 确定 4 的 值 ， 即 取 去 噪 图 像 质量 最 佳 


i 


EAE fe BR LEIA BAL) IA A 的 值 。 本 文 应 用 梯度 下 


降 法 求解 式 (10) 得 到 全 变 分 模型 去 噪 后 的 低频 子 图 像 ， 应 用 式 


(6) 


定 对 图 
HL, , €H AHT BE 


i= 


Er 


45 BF EE E 
小 波 逆 变换 得 到 去 噪 后 的 
体 步骤 如 下 : 
a) 小 波 变换 。 选 用 sym5 小 波 基 对 图 像 进行 小 波 变换 ， 确 
像 的 7 层 分 解 ， 得 到 低频 子 图 像 LL ， 水 平 高 频 子 图 像 
1 LH, 和 对 角 高 ; 


波 系数 分 别 进 


函数 处 理 后 的 高 频 子 图 像 ， 最 后 通过 


图 像 ， 方 法 实现 流程 如 图 


2 所 示 ， 


TEA E 


VE 
=H 


子 图 像 HH,， 


b) 全 变 分 处 理 。 将 低频 子 图 像 LL 用 TV 模型 进行 去 噪 处 


得 到 处 理 后 的 低频 子 图 像 LL 。 


c) 闵 值 化 处 理 。 对 三 个 方向 的 高 频 子 图 像 ， 按 式 (6) 对 其 小 


LH;. HH}. 


d) 图 像 重 构 。 对 应 用 


低频 部 分 LL, 和 高 频 部 分 


换 ， 


行 阔 值 化 处 理 ， 获 得 处 理 后 的 高 频 子 图 像 AIL,’ ~ 


TV 模型 和 小 波 阔 值 化 处 理 后 的 图 像 


HL; LH}. HH; 进行 小 波 逆 变 


且 进 行 重 构 ， 最 后 得 到 去 噪 后 的 图 像 结果 。 


定稿 K A, F: AP MAAS RRA ALR 


TT 3 ”实验 结果 与 分 析 
y 
小 波 变换 为 验证 本 文 提 出 方法 的 去 噪 性 能 ， 通 过 仿真 实验 对 本 文 提 
出 的 去 噪 方法 和 ROF 法 、 非 均值 滤波 (NLM) 去 噪 法 、 自 适应 小 
ae Sakae 水 平 高 频 系数 ET A E E RA KS 图 像 去 如 效果 进行 对 比 o 采用 Barbara、 Peppers 
和 Lena 三 种 不 同 细节 的 图 像 进行 实验 。 对 三 幅 图 像 加 入 均值 
Y AO, WAP RHE oO 分 别 为 5、20、25 和 30 的 加 性 高 斯 白 
TV 去 噪 自 适应 阐 值 去 噪 a ` Ro s= van ` y 
噪声 ， 为 确保 实验 可 行 性 ， 本 文选 用 sym5 作为 小 波 基 。 
3.1 去 噪 结果 对 比 
= 
图 像 重 构 为 观察 不 同 噪声 程度 的 去 噪 效果 , 这 里 将 加 入 标准 差 a 为 
y 15、20、25 和 30 的 高 斯 噪声 的 Barbara, Peppers 和 Lena 图 像 
分 别 用 ROF 法 、 非 均值 滤波 去 噪 法 、 自 适应 小 波 阔 值 去 噪 法 以 
图 2 ”去 噪 方法 实现 流程 图 及 本 文 提 出 的 去 噪 方法 进行 去 噪 的 结果 列 出 ， 如 表 1~3 所 示 。 


表 1 Barbara 图 像 不 同 噪声 程度 去 噪 结果 对 比 


ROF NLM Ai BY 改进 方法 
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20 


25 


30 
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30 


录用 定稿 K 弘 ， 等 : 基于 小 流 阅 值 和 全 变 从 模型 的 阁 像 去 噪 


表 3 Lena 图 像 不 同 噪声 程度 去 噪 结果 对 比 
oO ER 含 噪 图 像 ROF NLM EEM BRE 改进 方法 


15 


20 


25 


30 


为 了 更 好 地 比较 图 像 去 噪 前 后 结构 保持 情况 ， 本 文选 | 
Sobel 算 子 检测 用 ROF 法 、 非 均值 滤波 去 噪 法 、 自 适应 小 波 阔 
值 去 噪 法 以 及 本 文 提 出 的 去 噪 方法 分 别 对 加 入 均值 为 0、 标 准 
差 9 =30 的 高 斯 噪声 的 Barbara, Peppers 和 Lena 图 像 去 噪 结果 
的 边缘 ， 结 果 如 图 3~5 所 示 。 


ORREZ 


(b) 含 噪 图 像 


(c)ROF 法 (d)NLM 


(c)ROF 法 (d)NLM 


(e) 自 适应 小 波 闵 值 法 (HA SCT 


图 4 Peppers 去 噪 后 的 边缘 提取 图 像 

从 表 1~3 中 可 以 看 出 ， 随 着 噪声 方差 的 增 大 ， 采 用 ROF 
法 虽 能 抑制 一 定 的 噪声 ， 但 是 过 度 平 滑 了 图 像 中 的 细节 信息 ， 
如 Barbara 的 头巾 、 裤 子 和 桌布 、Lena 的 帽子 等 都 有 过 度 平滑 
适应 小 波 闵 值 法 (本 文 方法 WAR, ROF 进行 边缘 检测 时 ， 容 易 受 到 噪声 的 影响 ， 如 图 4(d) 
图 3 Barbara 去 噪 后 的 边缘 提取 图 像 中 Peppers 的 部 分 细节 丢失 。 并 且 从 表 3Lena 图 像 看 出 去 噪 后 


(e 


<~ 


录用 定稿 


KI 


不 
ARK 


后 


中 Lena 帽子 、_ Barbara 的 头 ! 
Barbara 的 头 部 和 Lena 帽子 
上 ， 本 文 方 法 去 噪 后 提取 的 边 
自 适 


ae 


J 


RIL 


BE 


背景 亮度 变 暗 ， 
缘 模 糊 ， 存 在 失真 现象 ， 如 
K 1 Barbara 的 头 部 和 手臂 。 采 
彻底 ， 且 同样 出 现 边 缘 模 尖 
景 以 及 Barbara 的 右 侧 桌布 均 有 模糊 现象 ， 


去 噪 效果 不 理 想 。 采 用 


K 


非 均值 ; 
3~5 Ja z 


Lena 


A} 


#1 


， 如 


Ai- 


图 像 平滑 自然， 


均值 滤波 法 和 


:具有 更 好 的 去 噪 效果 ， 


其 边缘 和 细节 


~3 中 Lena 头发 、 
用 本 文 方法 去 噪 
sth ae A 如 图 3~5 


适应 小 波 阔 值 法 


Dy 
7 


pa 


] 和 桌布 纹理 、Peppers 都 比较 清 


青 晰 ， 


的 轮廓 更 明显 ， 从 图 
缘 信息 更 完整 。 
应 阀 值 法 去 噪 结 果 的 对 
能 更 好 地 保护 


KI 


3~5(f) AT UA 
通过 与 ROF 法 、 
比 可 以 看 出 ， 


像 的 边缘 纹理 和 


从 本 


(c)ROF 法 


(d)NLM 


(6e) 自 适应 小 波 闵 值 法 HA SATE 
图 5 Lena 去 噪 后 的 边缘 提取 图 像 
3.2 评价 标准 
为 进一步 对 本 文 方法 进行 客观 的 评价 ， 这 里 选用 峰值 信 噪 


比 (peak signal noise ratio，PSNR)、 均 方 误差 (mean square error, 
MSE) 以 及 结构 相似 度 (structual similarity, SSIM) 


准 。 


体 公式 如 下 : 


255° MN 


m= 


[FG D-8, j) 


i=l j=l 


LFE j- gi, DP 


i 


EI 


F 价 标 


(12) 


(13) 


ik 5h, 


(2u,u, + c20 y +¢,) 


等 : 基于 小 波 阅 值 和 全 变 分 模型 的 图 像 去 噪 


ssm = (u +u,” 十 CI Yo,” +0,” +¢,) a 
其 中 : (OD RAARAR, 8G) 表示 去 噪 后 的 图 像 ，M 和 
NABI VAMT MAIR, Up, Up, OF, OF, We 
分 别 表示 图 像 / 、8 WE. TÆ. TH, ci ，c; 为 常数 。 
PSNR 表示 原 图 像 与 降 噪 图 像 噪声 水 平 的 差异 ， 能 够 反映 处 理 
后 图 像 逼 近 原 始 清晰 图 像 的 程度 。PSNR 值 越 大 去 噪 后 的 图 像 
在 平均 意义 上 越 接近 于 原始 图 像 。MSE 是 原始 图 像 和 去 噪 后 图 
像 之 差 平方 的 期 望 ，MSE 值 越 小 图 像 去 噪 效果 越 好 。SSIM 直 
接 将 两 幅 图 像 的 结构 相似 程度 作为 评价 标准 ，SSIM 值 越 大 ， 
说 明 去 噪 图 像 与 原始 图 像 越 相似 ， 去 噪 后 图 像 视觉 效果 越 好 。 
这 里 将 标准 差 设置 为 15、20、25 和 30， 计 算 Barbara, Peppers 
和 Lena 图 像 分 别 采取 自 适 应 小 波 阔 值 去 噪 法 、ROF 法 、 非 均 
值 滤波 和 本 文 改 进去 噪 方法 的 峰值 信 噪 比 、 均 方 误差 和 结构 相 
似 度 ， 观 察 噪 声 程度 对 去 噪 方法 的 影响 ， 其 结果 分 别 如 表 4~6 
所 示 。 
表 4 不 同方 法 的 峰值 信 噪 比 (dB) 对 比 
图 像 NBA ROF NLM 改进 方法 
79.34 78.00 73.90 80.70 
76.46 73.99 72.49 77.62 
Barbara 
74.51 71.20 71.82 75.71 
73.38 69.55 71.48 74.51 
75.22 77.44 79.05 81.88 
68.15 73.81 73.96 79.66 
Peppers 
65.92 71.12 70.74 78.03 
62.97 68.91 68.91 76.99 
73.96 76.50 79.11 82.89 
69.57 73.97 73.96 80.17 
Lena 
66.81 70.99 70.63 78.37 
63.78 69.57 68.91 77.26 
表 5 不 同方 法 的 均 方 误差 对 比 
图 像 o MBE ROF NLM 改进 方法 
15 0.74 2.51 2.21 0.93 
20 1.12 4.02 3.74 1.51 
Barbara 
25 1.77 6.30 4.12 1.71 
30 3.04 7.43 4.65 2.30 
15 0.32 1:21 0.58 0.17 
20 0.68 2.74 1.26 0.46 
Peppers 
25 1.07 5.03 1.53 1.02 
30 1.32 7.12 1.86 1.21 
15 2.20 1.01 0.95 0.62 
20 2.91 2.60 1.24 0.98 
Lena 
25 3.12 4.91 1.57 1.04 
30 4.23 6.42 1.92 1.30 


录用 定稿 
表 6 不 同方 法 的 SSIM 对 比 
图 像 oO E E BR) {EL ROF NLM 改进 方法 
15 0.809 0.824 0.802 0.834 
20 0.717 0.711 0.714 0.721 
Barbara 
25 0.698 0.673 0.683 0.703 
30 0.527 0.515 0.534 0.649 
15 0.867 0.811 0.825 0.874 
20 0.765 0.772 0.764 0.781 
Peppers 
25 0.641 0.643 0.635 0.652 
30 0.604 0.607 0.609 0.611 
15 0.813 0.809 0.819 0.822 
20 0.695 0.733 0.721 0.749 
Lena 
25 0.631 0.647 0.624 0.658 
30 0.607 0.625 0.617 0.631 
从 表 4 结果 可 以 看 出 ， 随 着 加 入 Ga 的 增 大 ， 四 种 方法 的 峰 
值 信 噪 比 都 逐渐 减 小 。 四 种 方法 中 ， 无 论 是 对 纹理 信息 丰富 的 


Barbara 和 Peppers 图 像 ， 


还 是 对 边缘 信息 比较 丰富 Lena 图 像 ， 


DA 


所 


E 


0 =15 时 得 到 的 SSIM 值 降 低 ] 
SSIM 值 平均 降低 了 0.234， 


用 本 文 方法 得 到 的 峰值 信 品 比 最 大 ， 均 方 误差 最 小 ， 如 表 5 
示 ， 说 明 本 文 方法 的 去 噪 效果 达到 了 最 好 。 

从 表 6 可 以 看 出 ， 当 标准 差 o =15 时 ， 四 种 方法 的 SSIM 
相差 不 大 ， 当 a =30 时 ， 本 文 方法 得 到 的 SSIM 值 相 比 于 
0.128, 而 其 他 三 种 方法 得 到 的 
随 着 O 的 增 大 本 文 方法 所 得 


SSM 值 最 大 ， 说 明 利用 本 文 方法 得 到 的 去 噪 图 像 视觉 效果 最 
好 。 
4 ARA 


Hy 


小 
分 
的 
本 


度 、 不 同方 向 上 的 含 噪 情况 确定 最 


针对 小 波 阔 值 方法 去 噪 不 彻底 且 导 致 图 像 边缘 模糊 的 问 
， 本 文 提 出 一 种 基于 自 适 应 小 波 阔 值 和 全 变 分 模型 相 结合 的 
像 去 噪 方法 ， 该 方法 应 用 小 波 阔 值 法 在 高 频 部 分 根据 不 同 尺 
RRM, ooh So 
波 统一 阔 值 的 不 足 ， 增 强 了 算法 的 自 适 应 性 。 同 时 利用 全 变 
法 对 低频 部 分 降 品 处 理 的 良好 效果 ， 更 多 地 保留 了 原始 图 像 
边缘 信息 。 对 含 不 同方 差 的 高 斯 噪声 图 像 的 实验 结果 表明 ， 
SC HE HH AY FTE MM BE HS A A RRE, OBOE 
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Dit 
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方 
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误差 和 较 高 的 峰值 信 噪 比 ， 理 细节 也 能 得 到 较 
的 保留 。 
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